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Abstract 

A scalar potential model (SPM) was developed from considerations of 
galaxy clusters and of redshift. The SPM is applied to Hi rotation curves 
(RCs) and RC asymmetry of spiral galaxies. The resulting model adds the 
force of a scalar potential of the host galaxy and of neighboring galaxies to 
the Newtonian rotation velocity equation. The RC is partitioned radially 
into regions. The form of the equation for each parameter of each region 
is the same with differing proportionality constants. Integer values of 
each equation are determined empirically for each galaxy. Among the 
sample galaxies, the global properties of galaxies of B band luminosity, of 
position, and of orientation determine the RC and RC asymmetry. The 
Source of the scalar field acts as a monopole at distances of a few kpc from 
the center of spiral galaxies. The scalar potential field causes Newtonian 
mechanics to considerably underestimate the mass in galaxies, which is 
the "missing mass problem". The SPM is consistent with RC and RC 
asymmetry observations of the sample spiral galaxies. 

Galaxies: kinematics and dynamics - Galaxies: fundamental parameters - 
Galaxies: interactions - Galaxies: spiral 



1 Introduction 



The discrepancy between the Newtonian estimated mass in spiral galaxies and 
the observation of star and gas kinematics is well established. The rotation curve 
(RC) relates the square of the tangential rotation velocity v (km s _1 ) of particles 
in orbit around spiral galaxies to their galactocentric radius R (kpc). The RC 
is measured along the major axis. Traditionally, the focus has been on account- 
ing for H I RCs that are flat in the outer region immediately beyond the knee 
(OR). However, observations also include rising RCs, declining RCs, an abrupt 
change in slope at the extreme outer region (EOR ) in many galaxies, and ro - 
tational asymmetry with non-relativist i c velo c ities. Battaner fc Florida (|2000|): 
iFerrarese fc Merrittl (l20o3):lGhez et ail (j2000h : ISoiue et all (|l999l) : ISofue fe Rubinl 
(1200 li) : and lTakamiva fc Sofuej (|2000f ) provide a summary of the knowledge of 
RCs. The RC is an excellent tool to evaluate galaxy models. 

The RC differs for different particles. For ex ample, the H I RC and the RCs 



of stars as shown by the Ha line for NGC 4321 (jSempere et all 119951 ) differ in 



the rapidly rising region before the knee (RR) and approach each other in the 
OR as shown in Fig. [TJ 

The stars near the center of galaxies have rotation velocities more than an 
order of magnitude larger than v at the knee of the Hi RC. The decline of the 
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Figure 1: Plots of the square of the rotation velocity v (10 km s ) versus 
galactocentric radius R (kpc) of the H i RC ( filled diamonds) and Ha RC (filled 



squares) for NGC4321 (jSempere et al.l . 1199 



RC (fil 
1995). 



stella r RCs from the center is Keplerian ( Ghez et all l2000l : iFerrarese fc Merritt , 
l2002h . The v continues to decline to the radius Ra (kpc) at which the rotation 
velocity va is at a minimum. The va is approximately an order of magnitude 
smaller than v at the knee of the Hi RC. 

The Hi RC rarely extends inward to near Ra- When drawing the Hi 
RC, a straight line is often drawn from the innermost measured point through 
the c enter of the galaxy. This practice is inconsistent with ste llar observa- 
tions ( Schodel et all 2002 ) and is incorrect ( Sofue fc Rubin . 2001 ). 

The evidence that the cluster environment affe cts RC shape appears contra- 
dictory ( Whitmore et all Il988t [Dale et al. . 2001 ). These studies concentrated 
on Ha RCs. The Hi RC usually extends farther outward than Ha RCs. There- 
fore, the Hi RCs are more sensitive to neighboring galaxy's effects. 

As larger samples of RCs have been measured, researchers attem pted to pro- 
duce a quantitative description of an empirically derived, average RC. lPersic et al 
(|l99d) claimed the existence of a universal RC (URC) for a large sample revealed 
by the feature that only luminosity of the h ost galaxy dictated th e RC. Indepen- 
dent work has confirmed much of the URC ( Catinella et al. , 20061 and references 
ther ein). Othe r researchers have discussed the i nadeq uacy of the average RC 
(e. g. ICourteaulll997t iGarrido et aDl2004t IWillickl [19991 



Roscoej (|2002l) used a dynamical partitioning process and found that the 
dynamics in the outer part of optical RCs are constrained to occupy one of 
four discrete dynamical classes. The classes are determined by the absolute 
magnitude, surface brightness, and a parameter for each optical RC that is an 
exponent of the radius at which v is measured. 

Galaxies with a small asym metry in the RC of a few km s _1 are well k nown 
dShane k. Bieger-Smitb , 19661 ) . Large RC asymmetries are less appreciated (jjosi 
20021 ). Because the observational data from each side of a galaxy RC are gen- 
erally averaged, only highly asymmetric cases are recognize d. RC asymmetry 
app ears t o be t he norm rather than the exception ( Jog , 2002t ). IWeinberg] (1995) 
and I Jod ( 1997 ) proposed the implied mass asy mmetry is d u e to an imposed 
lopsided potential caused by galaxy interaction. I Dale et al.l (|200lh found RC 
asymmetry of early type galaxies falls by a factor of two between the inner and 
outer regions of clusters. The formation, evolution, and long term maintenance 
of galactic, kinematical asymmetry remains a mystery. 

With the improved Hi measuring equipment in the 1990's, the Hi RC ob- 
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se rvations detected an a brupt change in the EOR in many galaxies (see figures 
in Catinella et al.l 20061) . This characteristic is ignored in current models. For 



example, iPizzella et al.l ([20051 ) considered this factor as a reason to exclude 
galaxies from the sample. 

Spectr a coming from HII regions are correlated systematically with R and 
little else (jBinnev fe Merrifieldl . ll.998l . pp. 516-522). The i nterstellar abun dances 
of metals in a disk galaxy decrease with increasing radius ( Tadrossl 12003 ). Also, 
the absolute B band magnitude Mb (mag.) of a galaxy is correlated with 
the metallicity obtained by extrapolating [O/H] within the disk to the galactic 
center. Low luminosity galaxies tend to be metal poorer than high luminosity 
galaxies. 

If most of the matter of a galaxy is in the bulge region, classical Newtonian 
mechanics predicts a Keplerian declining RC in the disk. The observation of 
rising RCs in the RR and of rising, flat, or declining RC's in the OR poses a 
perplexing problem. Postulates to model the RCs have included the existence 
of a large mass Mdm of non-baryonic, dark matter (DM) and dark energy; the 
variation of the gravitational acceleration R; the variation of the gravitational 
constant G with R; and the existence of an added term to the Newtonian grav- 
itation equation opposing gravitational attraction. 

Repulsive DM has been propo sed to provide a repulsive force that has the 
form of an ideal, relativistic gas (ICoodmanl . boOOh . In DM models the rising 
RCs of low surface brightness galaxies (LSBs) require larger relative amounts of 
DM than flat o r declining RCs. This contradiction to the successful Tully-Fisher 
relati on (TF) (jTullv fe Fisherl . Il977t ) is well established (|Sanders fc McGaughL 
20021 ). Another mystery of the DM paradigm is the apparent stro ng luminous 



to-da rk matter coupling that is inferred from the study of RCs (jPersic et al 
19961) . 



Currently, the most popular m odified gravity model is the Modified New- 
tonia n Dynamics (MOND) model ([Bottema et al.l . 120021 : ISanders fc McGaugb , 
2002L and references therein) . MOND suggests gravitational acceleration changes 
with galactic distance scales by a fixed parameter related to the mass-to-luminosity 
ratio with units of acceleration. MOND appears limited to the disk region of 
spiral galaxies. MOND requires the distance for NGC 2841 and NGC 3198 to be 



consi derably larger than the Cepheid-based distance D a (Mpc) ( Bottcm a et al 



2002). Appling MON D to cosmological scales is being investigated (jBekenstein 



2004; ISandersl .1 2006). MOND may r epresent a effective force law arising from 
a broader force law ( McGaugh . 19991 ). 

The possible variation of G or \G/G\ ~ H Q , where H a is the Hubble con- 
stant, were first proposed by Dirac. This was followed by the Brans-Dicke 
theory of gravity that postulated G behaves as the re c iproca l of a scalar field 



whose source is all matter in the universe (jNarlikarl . 120021 ). More recently 
iBrownstein fc Moffat ( 20061 . and references therein) have proposed a variation 
of gener al relativity in which G varies with R. T his model has been applied 
to RCs (jBrownstein fc Moffatl . 120061 : iMoffatl . 120051 ) and the postulated gravita- 
tional radiation of binary pulsars. This model has results similar to the MOND 
model. However, it us es a distance D (M pc) for NGC 2841 of 9.5 Mpc rather 



than D a =14.07 Mpc (jMacri et all . 12001). 



Alternate models of the RC posit the existence of another term in the Newto- 
nian equation resulting from a scalar field of a "fift h force" , a scalar field coupled 
to ma tter, or a static^ scalar field coupled to DM. IPadmanabhan fc Choudhurv 
( 2002L and references therein) explored the possibility of a scalar field acting 
as d ark matter. The usua l "fift h force " mod els posit a Yukawa like poten- 
tial (Bertolami & Paramosl 120051 ). iFavl (|2004l ) postulated minimally co upled 
and massive scalar fields are responsible for flat RCs, only. iMbelekl (2004) sug- 



1 Static in analogy to an electrostatic field caused by charged particles. 
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gcstcd a real scalar field, minimally coupled to gravity sourced by baryonic and 
dark matter may also reproduce flat and rising RCs with the assumption that 
the dark halo mass density dominates the baryonic mass density. The free pa- 
rameters are a turnover radius r D , the maximum rotation velocity, which is at 
r Q , and an integer. If a natural theoretical explanation could be found for the 
relation between r and the positive integer, the integer may be the only free pa- 
rameter. The postulate of an unknown "fifth force" is as likely as the postulate 
of unknown and undetected DM. 

The scalar potential model (SPM) was created to be consistent with the 
observation of the morphology-radius relation of galaxies in clusters, of the 
intragalactic medium of a cluster of galaxies, and of the flow of matter from 
spiral galaxies to elliptical galaxies ( Hodg e, 2006). The SPM was applied to 



rcdshift observations and found a qualitative explanation of discrete redshift 
observations. The SPM suggests the existence of a massless scalar potential 
p oc derived from a heat differential equation. Physically, the heat equation 
requires a flo'wH of energy from Sources to Sinks to form the potential field. That 
the Vp field acts like a "wind" was suggested. Several differences among galaxy 
types suggest that Sources of p are located in spiral galaxies and that Sinks of p 
are located in elliptical, lenticular, and irregular galaxies. The Source forming a 
galaxy leads to the proportionally of the Source strength and emitted radiation 
(luminosity). Therefore, the total mass of a galaxy is related to the luminosity 
of a galaxy. A cell structure of galaxy groups and clusters was proposed with 
Sinks at the center and Sources in the outer shell of the cells. The cell model is 
suppo r ted by the dat a and a nalysis of Aaronson et aD (19821) ; Ceccarelli et al. 
(|2005h : lHudson et all (|2004h ; iLilie et all (|l986f k and iRejkubal ((2004J). Because 



the distance between galaxies is larger than the diameter of a galaxy, the Sources 
were considered as point (monopole) Sources. 

The force F s (dyne) of the p field that acts on matter is 

F s =G s m s \7p, (1) 

where (1) the G s is a proportionality constant; (2) m s is the property of particles 
on which F s acts; and (3) p is the sum of the effects of all galaxies, 

Source A^ink 

where q and r\i are numbers representing the strength of the i th Source and I th 
Sink, respectively, K e and K v are proportionality constants, and (Mpc) and 
ri (Mpc) are the distance from a Source and a Sink, respectively, to the point 
at which p is calculated, e > 0, r\ < 0, and N source and are the number of 

So urces and Sink s, respectively, used in the calculation. 

Hod ge (2006) suggested the m s property of matter is the cross section of the 



particle; the e oc Mt e , where Mt e is the total mass of the Source galaxy; and the 
r\ cx Mt v , where Mt v is the total mass of the Sink galaxy. 

This Paper further constrains the SPM to describe another set of observa- 
tions inconsistent with present models. The SPM finds the RC in the OR may 
be rising, flat, or declining depending on the galaxy and its environment. The 
Source of the scalar field acts as a monopole at distances of a few kpc from 
the center of spiral galaxies. In addition to fitting the RR and OR parameters, 
the SPM is consistent with EOR and asymmetry observations. In section 2, 
the model is discussed and the SPM v 2 calculation equation is developed. The 
resulting model is used to calculate RC parameters in Section 3. The discussion 
and conclusion are in Section 4. 



2 A flow of energy is in contrast to a static field. An analogy is the flow of gas or fluid from 
Sources to Sinks. 
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2 Spiral galaxy model 

The coordinate system center was placed at the sample galaxy's kinematical 
center and was aligned to our line of sight. The SPM posits the simplified New- 
tonian calculation of v of a test particle is the first approximation of v in an 
isolated galaxy, that observed asymmetries and deviations from the simplified 
Newtonian calculation in spiral galaxies are the result of the influence of neigh- 
bor galaxies, and that the v 2 of a particle in orbit of a spiral galaxy is the sum 
of the effects of F s and the gravitational force F g . Therefore, 

^ = ^_ Gs I^_j|.a + |i?xa |, (3) 
R m L R 



wher e (1) the Newtonian simplifications of circular motion (R — 0) (jBinnev fc Merrifieldl . 



19981 p. 725) and a spherical mass density field are assumed, (2) the mass of 



the test particle is assumed to be constant over time, (3) M (M Q ) is the mass 
inside the sphere of radius R from Newton's spherical prop erty, (4) th e inertial 
mass m L equals gravitational mass m g of the test particle (|Willl . ll993l) . (5) the 



L term is due to the F s of the host galaxy, which also derives from a spheri- 
cally symmetric p field; (6) L = K e e = io -0 - 4Mb erg s _1 f or Source gal axies 
or L = K v r] = -2.7 x 1CT - 4M b erg s -1 for Sink galaxies dHoded . 12001 : (7) 
| indicates absolute value; and (8) a (km s~ 2 ) is the acceleration caused by 
neighboring galaxies, 

G s m s -± 

a Q = Vp, 4) 

where the number of galaxies exclude the host galaxy. Note that no assumption 
about the significance of a Q has been made. 

Because v is measured only along the major axis in the region under consid- 



eration (jBinnev &: Merrifieldl . Il998l p. 725) and if the Vp field is approximately 



uniform across a galaxy, Eq. becomes 

2 M G s m s L -* _ 

V = G— — h \K •a |.R Ina jor, (5) 

-^major -*£major 

where K (km kpc -1 ) is a constant vector and i? ma j r (kpc) is the galactocentric 
radius along the major axis. 

In Fig. [T] the H I RC at lower radius R„ (kpc) in the RR has two scalloped 
shapes that suggests spherically symmetric shells of matter. Also, the Ha RC 
rapidly increases, peaks, and then declines at the beginning of each shell. The 
Ha lines are generally formed in excited interstel lar gas. In the disk regio n 
of a galaxy, the gas is usually excited by hot stars (jBinnev fc Merrifieldl . [l998h . 



Because the m s /m L factor must be different for different matter types, each shell 
has a different metallicity star type. Because the Ha RC approaches the H I RC 
in the disk region such as plotted in Fig. [1] with hot, hydrogen burning stars, 
the m s /m L factor must be the same for Hi and hydrogen stars. This suggests 
the m s /m L factor varies by element type and acts on atoms at the largest. The 
metallicity - radius relation follows. 

The m s /m L ratio of stars is changing through changing elemental composi- 
tion by nucleosynthesis in addition to accretion and emission of matter. There- 
fore, the H I RC is preferred to trace the forces influencing a galaxy outside the 
bulge. Because only the Hi RC is considered in the calculations herein, the 
units used were G s m s /m L — 1 kpc km 2 s _1 erg -1 . 
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Table 1: Data for the select galaxies. 



Galaxy 


morphology" 


t b 


lc c 


rc" 




NGC 0224 


SA(s)b 


3 


2 


F 


0.79 


NGC 0300 


SA(s)d 


6.0 


6 


R 


2.00 


NGC 0598 


SA(s)cd 


6 


4 


R 


0.84 


NGC 0925 


SAB(s)d HII 


7 


4 


D 


9.16 


NGC 2403 


SAB(s)cd 


6 


5 


F 


3.22 


NGC 2841 


SA(r)b;LINER Sy 


3 


1 


F 


14.07^ 


NGC 3031 


SA(s)ab;LINER Syl.8 


2.4 


2 


D 


3.63~ 


NGC 3198 


SB(rs)c 


5.2 


3 


F 


13.80 


NGC 3319 


SB(rs)cd;HII 


5.9 


3.3 


R 


13.30 


NGC 4258 


SAB(s)bc;LINER Syl.9 


4 


3.9 


D 


7.98 


NGC 4321 


SAB(s)bc;LINER HII 


4 


1 


D 


15.21 


NGC 4414 


SA(rs)c? LINER 


5.1 


3.6 


D 


17.70 


NGC 4535 


SAB(s)c 


5 


1.9 


D 


15.78 


NGC 4548 


SBb(rs);LINER Sy 


3.1 


2 


D 


16.22 


NGC 5457 


SAB(rs)cd 


5.9 


1 


D 


6.70 


NGC 7331 


SA(s)b;LINER 


3.9 


2 


F 


14.72 



"Galaxy morphological from the NED database. 

'Galaxy morphological type code from the LEDA database. 

c Galaxy luminosity class code from the LEDA database. 

^Galaxy's Hi RC type according to slope in the OR. R is rising, F i s flat, and D is declining. 
e The distance £> a (Mpc) to the gala xy from lFreedman et alj <|200lT > unless otherwise noted. 
^The D a is from lMacri et al.l [120011) . 



3 Results 
3.1 Sample 

The elliptical, lenticular, irregula r , and spiral galaxies used in the calculations 
were the same as used in iHoded (120061, . That is, they were selected from the 



NED databas<§ The selection criteria were that the heliocentric redshift z mn 
be less than 0.03 and that the object be classified as a galaxy. The parameters 
obtained from the NED database included the name, equatorial longitude £ion 
(degrees) for J2000.0, equatorial latitude Ei at (degrees) for J2000.0, morphology, 
the B-band apparent magnitude nib (mag.), and the extinction E x t (mag.) as 
defined by NED. The galactocentric redshift z was calculated from z mn . 

The 21-cm line width W20 (km s _1 ) at 20 percent of the peak, the inclination 
i n (arcdegrees), the morphological type code "t", the luminosity class code "lc", 
and the position angle P a (arcdegrees) for galaxies were obtained from the LEDA 
database] when such data existed. 

The selection criteria for the "select galaxies" we re that a H I RC was av ail- 
able in the l i teratu re and that £> a was available from Freedman et al. ( 200 ll ) or 



Macri et al.l (|200ll ). Data for the 16 select galaxies are shown in Tableland 
their H 1 RCs arc plotted in Fig. [2] The L for the select galaxies was calculated 
using D a , TOb, and E xt . The selection criteria for the other sample galaxies were 
that a H 1 RC was available in the literature and that the L could be calculated 
using the TF method with the constants developed in [Hodge (|2006h . 



This select galaxy sample has LSB, medium surface brightness (MSB), and 
high surface brightness (HSB) galaxies; has LINER, Sy, HII, and less active 
galaxies; has galaxies that have excellent and poor agreement between the dis- 
tance D t f (Mpc) calculated using the TF relation and D a ; has a D a range of 
from 0.79 Mpc to 17.70 Mpc; has field and cluster galaxies; and has galaxies 
with rising, flat, and declining RCs. 



3 The Ned database is available at http://nedwww.ipac.caltech.edu The data were ob- 
tained from NED on 5 May 2004. 

4 The LEDA database is available at http:/ /leda. univ-lyon.fr The data were obtained from 
LEDA on 5 May 2004. 
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Figure 2: Plots of the square of the H I rotation velocity v 2 (10 3 km 2 s -2 ) versus 
galactocentric radius i? m ajor (kpc) along the major axis. The straight lines 
mark the application of the derived equations to the RCs of the select galaxies. 
The application of the derived equations to NGC 0224, NGC 0300, and NGC 
0598 were omitted because these galaxies lacked a |i^»a | val ue. The references 
for th e RCs are noted in bracke ts and a re as follows: [1 Gottesman et al 



(Il966h. f2lCarignan fc Freeman! (Il985h. r3lCorbelli fe Salucci 



[4; 

flOlvan Albada fc Shand (11975D. flllSempere et al 



Krumm fc Salpeteij (Il979l ). [5u3egeman et alj (Il99ll). f6lRotd 



van Albada et all dl985TK [8lBosmal (Il98ll). f9lMoore fe Gottesmanl 



2000), 



1975) 



1998), 



19951). [12lBraine et al 



(Il993l). fl3lChincarini fc de Souzal (I1985D. \U 
[15 ]Huchtmeierl (|l975l ). and [16 ]Mannheim fc Kmetkol ( 



Vollmer et all (|l999f ). 
19961) . 
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3.2 First approximation 

When modeling an RC, the knee or turnover of an H I RC is often drawn smooth 
and relatively rounded. As depicted in Fig.[2]most of RCs of the selected galaxies 
have an abrupt change of slope at the turnover. 

The end of the RR was found by fitting a least squares straight line to the 
first three data points, if available, inward and immediately before the knee of 
the RC with a 0.99 or higher correlation coefficient. NGC 2841 lacked data 
points to identify the RR. NGC 4414 had only two data points before the knee 
and a significant decline in v 2 to indicate the RR. The linear relation suggests 
a cylindrical distribution of matter in the RR near the knee. 

At i?A, va ~ that implies the effective mass w inside a sphere with a 
i?A radius. That is, from Eq. (J5]) Ma ~ G s m s L/Gm L: where Ma is the mass 
within i?A at which dv 2 /dR = and at which the RR begins. 

The mass in an elemental volume dR in the RR near the turnover was 



modeled as a cylinder shell of small height H rr (|Binnev fc Merrifieldl . 119981 p. 
724) and with density D rr at i? rr . The mass M rr within a sphere with a i? rr 
radius is 

M rr = (D„H„2k) R 2 t + M A + M c , (6) 

where M c is a small correction term, the mass in the thin cylindrical RR shell ~ 
the mass in the spherical shell with the same radial parameters, and the Ma is 
distributed in a much thicker volume, which may be spherical. The components 
of M c are — (D lr H lr 2ir) R\, the small amount of matter indicated by va ^ 0, 
and the variation of the mass in the RR from a cylindrical distribution. 
Inserting Eq. (J6]) into Eq. ([5|) for the RR yields 

■j _ / m s. 



< = ^-*-GD„H„2w j Rrr 

G(M A + M c ) - m s G s L/m L 

-\K*CL \R, t , (7) 

where v rl (km s _1 ) is the rotation velocity of Hi in the RR at R IT . 

By averaging from side-to-side of the galaxy, the effect of the \K»a Q \ term was 
reduced and was considered insignificant for the first approximation compared 
to the D TT term near the knee. The F g and the F s offset each other at Ra to pro- 
duce a minimum v that then rises in the RR. Because «a ~ and the RCs near 
the knee are well fit by the straight lines, the [G(Ma t + M c ) — m s G s Lm~ 1 ]/ R„ 
term is small and nearly constant in the RR compared to the D TT term near the 
knee of the RC. Equation ([7|) suggest v 2 T depends on galaxy parameters D IT and 
H rr . Therefore, the first approximation is 

S TT - h 7 rr , (8) 



10 3 km 2 s- 2 kpc 

where S TI and / rr are the slope and intercept of the linear v 2 T — R lr relationship, 
respectively. 

As R rT increases, m s /m L of the predominant mass increases and the pre- 
dominant particles have decreasing density. Decreasing the amount of mass in 
a given cylindrical shell causes the end of the RR and the beginning of the 
transition region (TR) between the RR and OR. Because e cx M te and e cx Ma, 
the mass M rrmax within the maximum RR radius i? rrm ax (kps) oc e cx L of 
Eq. ([5|). The rotation velocity t) rrmax (km s _1 ) of the galaxy at i? rrmax can be 
used to compare parameters among galaxies. Because this change is caused by 
the change in elemental types that depends on the e and, hence, L of a galaxy, 
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the end of the RR is occurring under similar conditions in all galaxies. Thus, 
the f 2 rmax an d -Rrrmax are comparable parameters among galaxies. Posit, as a 
first approximation, 



v 2 

u rrmax 



103km 2 s-2 108 erg s-i +/a ' (9) 

where S a and I a are the slope and intercept of the linear relation, respectively. 

When the error in calculating distance is a small factor such as when using 
D a or D t [ , more galaxies were sampled. The dat a of v 2 rmav versus L of 1 5 select 
galaxi es from Fig. [21 and of 80 oth e r galaxies from Begeman et al. (Il991 ^ Broeils 



(Il992h; [Garcfa-Ruiz et alj (l2002h; iGuhathakurta et alj d 19881): iKornreich et al 



(2000. i2001h:lLiszt fc Dickevi (119951) : iMannheim fc Kmetkol (|1996l ): iRubin et al 
19851 ): and lSwaters et al.l ( 19991 ). for which the RR was identified, are plotted in 



Fig. [3j The other sample galaxies are listed in Table [2] denoted by an integer in 
the a\ column. Some of the galaxies had only two data points and a significant 
decline in v 2 to indicate the end of the RR. The distribution of the data points 
in Fig. [3] suggested a grouping of the galaxies such that the v 2 lmax versus L 
relations are linear for each group. This is reflected in the plotted lines in 
Fig. [3] Call each group a "classification". For the calculation, the (L,v 2 tmax ) = 
(0,0) point was included in all classifications. The relationship of the slopes of 
the lines is 

where a\ is an integer, S cvw d = 1.3±0.2 and / e rvd = 0.31 ±0.01 were obtained by 
a least squares fit to produce the lowest uncertainty a c . The a e was calculated 
as the standard deviation of the relative difference (^ 2 rmax /v 2 rmax ) between 
the measured parameter value and the calculated parameter value divided by 
the calculated parameter value for the sample galaxies. Because the goal is to 
examine the slope of the lines and because the (L,v 2 max ) — (0,0) point was 
included in all lines, minimization of er c chooses the calculated line with the 
closest slope to the slope of the line from (0,0) to the data point. A \ 2 type of 
minimization would choose a different line for some of the sample galaxies. 
Therefore, 

S a = K ai BH , (11) 

where I a = and K ai = 1.3 ± 0.2 and B ai = 2.06 ± 0.07 at 1 a. The values of 
K ai and B &1 were chosen as the maximum B &1 value that yielded a correlation 
coefficient for each of the lines from the origin to the data points (see Table [3] 
and Fig. [3]) greater than 0.90, that yielded a minimum <r e , and that rejected the 
null hypothesis of the test in Appendix A with a confidence greater than 0.95. 
Combining Eqs. ([9]) and (fTT|) yields 

„2 J 

103km 2 s- =iraiS -I^ri^ ±a °' (12) 

where a c = 21%. Note tr c = 16% for only the select galaxies. The large, filled cir- 
cle in Fig.[3]denotes the data point for NGC 5448 (Sv 2 lmeiy ./v 2 ] . max = —0.35). The 
large, filled square in Fig.[3]denotes the data point for NGC 3031 (£v 2 rmax / w 2 rmax = 0.25) . 

Tables [2] and [4] lists the a\ values for the sample galaxies. Table [5] lists the 
minimum correlation coefficient C cm in of the lines, the constants i4T x , and the 
exponential bases B x , where the subscript "x" denotes the generic term. 

The average difference between L and the luminosity calculated using D t { 
is 0.38L for the select galaxies. The relative difference in i> 2 rmax included the 
measurement uncertainty and the uncertainty that the RR may extend farther 
than the measured point such as seen for NGC 4258. The relative differences 
SL/L = 0.3 and <5v 2 rmax /v 2 rmax = 0.1 are shown as error bars in Fig. [3] for the 
select galaxies. 
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Table 2: Interger values for the non-select galaxies in the sample. 



Galaxy 


ai 


«2 


h 


h 


h 


h 


Galaxy 


ai 


a 2 


A 


/2 


h 


J2 


IC 0467 


3 


9 










NGC 


4206 


2 


7 










IC 2233 


4 


6 


4 





1 


5 


NGC 


4216 


2 


8 










NGC 0701 


3 


8 










NGC 


4222 


5 


10 










NGC 0753 


4 


8 










NGC 


4237 


4 


6 










NGC 0801 


5 


8 










NGC 


4359 


4 


3 


3 


1 


1 


5 


NGC 0991 


2 


3 


2 


-3 


1 


3 


NGC 


4378 


3 


11 










NGC 1024 


4 


8 










NGC 


4388 


5 


5 










NGC 1035 


4 


6 










NGC 


4395 


6 


3 


6 


-1 


-1 


5 


NGC 1042 


3 


4 


3 


-2 


3 


3 


NGC 


4448 


4 


7 










NGC 1085 


5 


6 










NGC 


4605 


4 


10 










NGC 1087 


4 


8 










NGC 


4647 


3 


7 










NGC 1169 










4 


5 


NGC 


4654 


4 


1 










NGC 1325 


2 


9 










NGC 


4682 


3 


7 










NGC 1353 


5 


6 










NGC 


4689 


5 


9 










NGC 1357 


4 


9 










NGC 


4698 


4 


6 










NGC 1417 


5 


9 










NGC 


4772 










1 


3 


NGC 1421 


3 


4 










NGC 


4845 


3 


7 










NGC 1515 


5 


4 










NGC 


4866 


4 


9 


4 


7 


3 


4 


NGC 1560 


4 


10 


5 


2 






NGC 


5023 


4 


8 


4 


-1 


1 


4 


NGC 1620 


3 


5 










NGC 


5107 













4 


NGC 2715 


2 


7 










NGC 


5229 


4 


6 


3 


4 


2 


4 


NGC 2742 


3 


5 










NGC 


5297 


4 


6 








2 


4 


NGC 2770 










3 


4 


NGC 


5301 


4 


8 


3 


5 


3 


3 


NGC 2844 


3 


13 










NGC 


5377 










2 


3 


NGC 2903 


5 


11 


3 


5 






NGC 


5448 


5 


7 


3 


3 


1 


2 


NGC 2998 


3 


2 










NGC 


5474 


1 


12 


1 


-2 





4 


NGC 3054 


4 


11 










NGC 


6503 


5 


9 


5 


4 






NGC 3067 


5 


15 










NGC 


6814 


7 


6 


8 


-3 






NGC 3109 


5 


5 


5 


2 






NGC 


7171 


3 


6 










NGC 3118 


4 


4 


4 


1 


2 


3 


NGC 


7217 


4 


8 










NGC 3200 


4 


15 










NGC 


7537 


4 


-1 










NGC 3432 


4 


10 


4 


6 


4 


5 


NGC 


7541 


4 













NGC 3495 


2 


5 










UGC 


01281 


4 


5 


4 


-4 





3 


NGC 3593 


5 


11 










UGC 


02259 


5 


8 


5 


1 






NGC 3600 


4 


5 


4 


6 


4 


5 


UGC 


03137 


5 


9 


3 


4 


3 


5 


NGC 3626 










4 


6 


UGC 


03685 


4 


7 


5 


-4 





6 


NGC 3672 


3 


8 










UGC 


05459 


5 


6 


4 


2 


2 


6 


NGC 3900 










1 


3 


UGC 


07089 


4 





4 


-2 





4 


NGC 4010 


3 


8 


4 


3 


2 


5 


UGC 


07125 


4 


4 


4 


-1 


1 


3 


NGC 4051 


3 


7 


1 


3 






UGC 


07321 


2 


7 


4 


8 


6 


4 


NGC 4062 


4 


5 










UGC 


07774 


4 


8 


4 


1 


2 


5 


NGC 4096 


4 


7 


3 


-2 


2 


4 


UGC 


08246 


3 


2 


3 


-2 





5 


NGC 4138 










2 


5 


UGC 


09242 













23 


NGC 4144 


3 


6 


4 


1 


2 


6 


UGC 


10205 


1 


10 










NGC 4178 


4 


6 
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100 




L (10 8 ergs" 1 ) 

Figure 3: Plots of square of the rotation velocity w„ max (10 3 km 2 s~ 2 ) at the 
maximum extent of the RR versus B band luminosity L (10 8 ergs -1 ) for the 95 
sample galaxies. The 15 select galaxies have error bars that show the uncertainty 
range in each section of the plot. The error bars for the remaining galaxies are 
omitted for clarity. The straight lines mark the lines whose characteristics are 
listed in Table [31 The large, filled circle denotes the data point for NGC 5448. 
The large, filled square denotes the data point for NGC 3031. 



Table 3: Data for the integer values of a\ of Eq. (fT2|) shown as the plotted lines 
in Fig. [3] 



„ a 

ai_ 


No." 


Corr. c 


F Test 


c> de 








1 


2 


1.00 


1.00 


2.482±0.004 


0.006±0.006 


2 


7 


0.97 


0.93 


6.6 


±0.6 


0.1 


±0.7 


3 


20 


0.99 


0.96 


12.8 


±0.5 


-0.8 


±0.8 


1 


.39 


0.98 


0.92 


23.2 


±0.7 


0.7 


±0.7 


5 


24 


0.97 


0.90 


42 


±2 


1 


±2 


6 


2 


1.00 


1.00 


84.3 


±0.3 


-0.2 


±0.2 


7 


1 


1.00 


1.00 


261 










"The integer denoting the exponential value in Eq. 11121 . 

'The number of data points for each ai value. 

c The correlation coefficient rounded to two decimal places. 

d The least squares fit slope of the n^ rmax - L lines. 

e The lines were calculated using (L,i)^ rmax ) = (0,0) as a data point. 

^The least squares fit intercept of the fj? rmax - L lines. 
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Table 4: First approximation integer values for the select galaxies. 



Galaxy 


at 


h 


Cl 


di 


ei 


fi 


5i 


ft! 


n 


h 


ki 


Ji 


NGC 0224 


5 


5 


5 


3 


4 


4 


3 


3 


4 






4 


NGC 0598 


4 


G 


3 


G 


G 


5 


5 


3 


G 






G 


NGC 3031 


5 


5 


5 


5 


4 


3 


3 


2 


4 


3 


5 


7 


NGC 0300 


5 


7 


3 


G 


7 


5 


5 


2 


7 






7 


NGC 2403 


5 


G 


4 


G 


G 


4 


4 


2 


G 


1 


4 


G 


NGC 5457 


4 


5 


5 


1 


2 


1 


4 


3 


1 


4 


5 


4 


NGC 4258 


3 


2 


1 


4 


5 


2 


2 


2 


3 


2 


1 


2 


NGC 0925 


5 


G 


4 


4 


2 


2 


4 


2 


4 


2 


3 


G 


NGC 2841 




1 








4 


4 


4 


3 


2 


5 


1 


NGC 3198 


4 


G 


5 


2 


1 


2 


4 


3 


3 


2 


G 


5 


NGC 4414 


5 


4 


4 


4 


4 


2 


3 


2 


3 


1 


2 


6 


NGC 3319 


3 


G 


4 


2 


3 


3 


4 


2 


5 


1 


5 


4 


NGC 7331 


5 


3 


3 


5 


5 


3 


3 


3 


4 


3 


3 


4 


NGC 4535 


4 


4 


3 


3 


4 


2 


1 


1 


5 


3 


4 


3 


NGC 4321 


4 


4 


4 


2 


2 


2 


1 


2 


4 


3 


2 


3 


NGC 4548 


6 


1 


2 


10 


8 


4 


2 


2 


G 


9 


1G 


7 



Other parametric relationships to L in the RR were calculated using the 
same procedure that was used to evaluate the ^ rmax - L relation. Because these 
equations involved -Rmnax, the calculation considered only the 15 select galaxies. 
The resulting equations are: 

^± = K ± u% ; (13) 

kpc 1 bl 10 s ergs- 1 ' v ' 



R v 2 

^rrmax / jjrrmax _ ts n di 



Im "°™ = K C1 B c c ] ± 15% ; (14) 

10 3 kpckm 2 s- 2 1 Cl lQSergs" 1 v ' 



K di K] T7T«— "I ± H% ; and (15) 



lO^kpc^kniV 2 1 dl 10 8 ergs 

10»kpc-\^-3 " K " B " TO^FT ± ** ; (16, 

where first approximation integer values are listed in Table 2] and the C cm i n , 
K K , and B x are listed in Table [5] 

Equations (fr2 |l - (fT5|) suggest at least two of the four integers ai, bi, ci, and 
di are functions of the other two. 

The end of the EOR was considered to be at the largest radius -Reormax (kpc) 
along the major axis that Hi is in orbit around the galaxy. Several galaxies 
have a H I measured point on one side of the galaxy at a greater i? ma jor than the 
other side. Because averaging the v from each side of the galaxy was required 
to reduce the \K • a Q \ effect that causes non-circular rotation, the outermost 
measured point used in the calculations was at the i? C ormax with data points on 
both sides of the galaxy. At R > -Reormax particles are no longer in orbit. For 
particles other than Hi, the m a /m b is less. Therefore, their maximum radius 
in a galaxy is less. N GC 4321 has Hi exte nding farther than a radius of four 



arcminutes. However. iKnapen et alJ (|1993f ) concluded the Hi motion beyond 



four arcminutes is not disk rotation. Therefore, -Reormax was considered to be 
at four arcminutes for NGC 4321. 

Because e tx M te and e oc Ma, Af oormax oc e tx L, where M eormax is the 
mass within i? cor max of Eq. ([5]). The rotation velocity iWmax of the galaxy at 
a radius of -Reormax can be used to compare parameters among galaxies. 
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Table 5: Values of the minimum correlation coefficients C cm i n , constants K, 
and exponent bases B x for the first approximation equations. 



Integer™ 


r ■ b 


K 




B*± 


Ol 


0.97 


1.3 


± 0.2 


2.06±0.07 


6i 


0.86 


0.34 


± 0.04 


1.88±0.05 


Cl 


0.96 


19 


± 1 


1.89±0.05 


tii 


0.99 


1.0 


± 0.1 


1.57±0.06 


ei 


0.97 


0.49 


± 0.06 


1.71±0.04 


h 


0.96 


8.6 


± 0.7 


1.57±0.04 


91 


0.93 


4.4 


± 0.4 


1.77±0.06 


hi 


0.97 


150 


±10 


2.5 ±0.1 


h 


0.95 


0.09 


± 0.02 


1.9 ±0.2 


h 


0.96 


0.075± 0.009 


1.95±0.02 



"The integer of the exponential value of parametric relationship that denotes the applicable 
equation. 

'The minimum correlation coefficient rounded to two decimal places of the lines of the 
parametric relationship. 

c The constant of proportionality at la of the parametric relationship. 
d The exponent base at la of the parametric relationship. 



Parameter relations to L in the EOR were calculated using the same proce- 



dure used to evaluate the 



L relation. 



The data of vl ve r sus L of 16 se l ect galaxies and of 34 ot h er galaxies 
fromlBegeirian et alj 3l99lh:lBroeild dl992h:lGarcia-Ruiz et all (l2002l) : lKornreich et al 
(|200dl200ll ): lLiszt fc Dickevl (|l995l ): lMannheim k Kmetkol (| 19961 ): and liwaters et &\ 



19991 ) for a total of 50 sample galaxies are plotted in Fig.|U The 34 other galax- 
ies are listed in Table [2] denoted by an integer in the /i column. The uncertainty 
in determining the value of the outermost point of the EOR is greater than de- 
termining the end of the RR. Because the determination of the outermost value 
of the parameters requires instrumentation sensitive enough to measure the val- 
ues, data prior to 1990 of 46 other galaxies used in the i>r rmax -L relation were 
omitted. The equations involving i? eormax used only the 16 select galaxies. 
The resulting equations are: 



10 3 km 2 s" 2 



10 8 ergs" 



— ± 13% , 



(17) 



where the large, filled circle in Fig. [4] denotes the data point for NGC 5448 
(^eormax/^oormax = — 0.20) and the large, filled square denotes the data point 



for NGC 3031 (6v; 



R 



2 

cormax 



= -O.i 



R, 



cormax 



kpc 



L 



1 10 8 ergs 



— ± 16% ; 



cormax "cormax 

10 3 kpckm s -2 

cormax 

corm 



L 



10 8 ergs _1 
L 



± 19% : and 



± 15% ; 



(18) 



(19) 



(20) 



H^kpc^WV 2 " ll 11 lOSergs- 1 

where first approximation integer values are listed in Tables [5] and 2] and the 
Ccmin) K x , and B x are listed in Tabled 

The uncertainty in calculating |if«a | depends on the accuracy of measuring 
the distance D (Mpc) to neighbor galaxies that are without a D„ and with- 
out a D t f measurement. The cell structure of galaxy clusters (lAaronson et al. , 
1982HCeccarelli et all l2005HHodgel . [2001 iHudson et all 12004 iLilie et all 11986 : 
Reikubal 12004) suggests a systematic error in z relative to D of a galaxy. There- 
fore, using z and the Hubble Law to determine D introduces a large error. 
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Figure 4: Plots of the square of the H I rotation velocity i>e 0rmax (10 3 km 2 s~ 2 ) at 
extreme outer region versus B band luminosity L (10 8 ergs -1 ) for the 50 sample 
galaxies. The 16 select galaxies have error bars that show the uncertainty level 
in each section of the plot. The error bars for the remaining galaxies are omitted 
for clarity. The large, filled circle denotes the data point for NGC 5448. The 
large, filled square denotes the data point for NGC 3031. 



However, the cell model also suggests neighboring galaxies have a similar D/z 
ratio in the first approximation. For those few galaxies near our line of sight to 



ratio m the nrst approximation, .tor those tew galaxies near our line ot sigiit to 
the sample galaxy and near the sample galaxv. lHodgel (2006) found a z change 



caused by the light passing close to galaxies. For a single galaxy this effect is 
small and was ignored. Therefore, 

A - —D p , (21) 

Zp 

where z p and D p (Mpc) are the redshift and distance of the sample galaxy, 
respectively, and Zi and Di (Mpc) are the redshift and calculated distance of 
the i th neighbor galaxy, respectively. The D p — D a for the select galaxies and 
D p = D t f for the other sample galaxies. The L of the i th neighbor galaxy was 
calculated using Di, mb, and E xt - 

Because larger \zi — z p \ implies larger error in Di, the N somccs and iV s i n ks of 
Eq. ([4| was limited to the number N of galaxies with the largest influence on 
|Vp| of the sample galaxy, where the Vp of Eq. (j4|) is evaluated at the center 
of the sample galaxy. The N = 7 was chosen because it produced the highest 
correlation in the calculations. Also, 6 < N < 9 produced acceptable correlation 
in the calculations. 

The calculation of \K • a D | also requires knowledge of the orientation of the 
sample galaxy. The direction of the sample galaxy's polar unit vector e po iar was 
defined as n orthward f rom t h e center of the sample ga laxy along the polar axis o f 
the galaxy. iHu et alj (|2006fk IPereira fc Kuhnl (|2005l ); and lTruiillo et all (|2006l ) 



found an alignment of the polar axis of neighboring spiral galaxies. Therefore, 
the orientation of e po iar with the higher |e po iar • o. a \ was chosen for the calcula- 
tions. 

The major axis unit vector e ma j r was defined as eastward from the center 
of the sample galaxy along the major axis of the galaxy. The minor axis unit 

Vector Grninor = ^major X 6p lar- 

At equal i? ma jor > -Rnmax on opposite sides of the galaxy, the L terms of 
Eq. (|5|) are equal and the M terms are nearly equal. Define asymmetry A sym , 

A sym = (v 2 h -vf)\ R , (22) 
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70 




IK-^I (10 8 kpc 1 km 2 s' 2 ) 4 

Figure 5: Plots of maximum asymmetry A symax (10 3 km 2 s~ 2 ) versus \K * a Q \ 
(10 3 kpc -1 km 2 s~ 2 ) for the 50 sample galaxies. The 13 select galaxies have error 
bars that show the uncertainty level in each section of the plot. . The error 
bars for the remaining galaxies are omitted for clarity. The large, filled circle 
denotes the data point for NGC 5448. The large, filled square denotes the data 
point for NGC 3031. 



where i>h and v\ are the larger v and smaller v at the same i? ma j r Qr) on 
opposite sides of the galaxy, respectively. 

The A sym is a function of -R ma jor- The A sym and \K • a \ are sensitive 
to both radial (non-circular) motion of particles and to variation of v due to 
torque. For the first approximation the maximum asymmetry ^4 symax in the OR 
depends predominately on \K •a \. Because Ajymax is minimally dependant on 
R, ^symax is comparable among galaxies. 

NGC 0224, NGC 0300, and NGC 0598 were omitted from the sample be- 
cause their z < and their neighbor galaxies may have positive or negative 
z. Therefore, evaluating the distance of neighbor galaxies of these galaxies was 
considered unreliable. Table [5] lists the \K • a a \, A symax , and A symax references 
for the remaining 13 select g alaxies. The sample consi st ed of 13 select g a laxies 
and o f 36 other galaxies from Garcfa-Ruiz et a l. ( 2002) ; iKornreich et al. ( 2000l 



2001 ): and Swaters et al.1 ( 1999h for a total of 49 sample galaxies. The 36 other 



galaxies are listed in Table [2] denoted by an integer in the ji column. The plot 
of j4 S ymax versus \K • a Q \ is shown in Fig. 

The same procedure that was used to evaluate the u 2 rmax - L relation was 
used to evaluate the A symax - \K • a \ relation. The result is 

^ symax - K h B\' \ K ' *o\ ± 28% ^ (23) 



10 3 km 2 s- 2 Jl jl lOSkpc-^m^^s^ 2 

where first approximation integer values are listed in Tables [2] and QJ C cm i n = 
0.98, K h = (1.000 ± 0.001), B h = 2.94 ± 0.09, and 

K = (3.60 x 10~ 6 ) e po i ar 
-(6.81 x 10- 6 )e major 

-(1.18 x 10- 5 )e minoT kmkp C - 1 (24) 

at la. The value of K was chosen such that K- n w 1.00. 

The large, filled circle in Fig. O denotes the data point for NGC 5448 
(<5A symax /A symax = —0.08). The large, filled square denotes the data point for 
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Table 6: Asymmetry and \K • a Q \ data for the select galaxies. 



Galaxy 


A a 


\K • aJ b 


Ref/. 


NGC 0925 


12.5 


1.53 


[11 

L J 


NGC 2403 


2.0 


0.99 


[2] 


NGC 2841 


4.0 


0.52 


[3] [4] 


NGC 3031 


7.1 


0.35 


[4] 


NGC 3198 


2.9 


0.27 


[5] [6] 


NGC 3319 


1.4 


0.44 


[7] 


NGC 4258 


14.0 


2.0 


[8] 


NGC 4321 


64.7 


2.6 


[9] 


NGC 4414 


6.5 


2.2 


[10] 


NGC 4535 


16.7 


0.58 


[11] 


NGC 4548 


27.9 


0.0018 


[12] 


NGC 5457 


34.8 


0.37 


[13] 


NGC 7331 


17.6 


0.71 


[6] 



a The units are 10 3 km 2 
The units are 10 3 kpc 
References: [lj Krum m 



2 . The uncertainty is ±10%. 
1 km 2 s- 2 . T he un cert ainty is ±30%. 



Salpeterl ll"l979h. [2lBegeman et ail (ll99lT). [31 iFillmore et al 



1986ft . [4 lRotsl dl975fl. f5lvan Albada et al] <198^7T6iBosmal il98lf >. f7 lMoore fc Gottesman 
19981). f8lvan Albada fc Shand J1975I1. [9tSempere et al I (|l995l 1. flOlBraine et al] [119931 1. 



fll lChincarini & de Souzal 1 119851) . [12 lVollmer et al.l Jl99gf) . and [13 lBosma et al.l lll98lf ) 



NGC 3031 (<L4 sy maxAWmax = -0.21). For the 13 select galaxies, a c = 17% 
[excluding the (A symax , \K • a \) = (0, 0) point]. 

Because the net effect of F s and F g from the host galaxy is weak at the outer 
edge of the galaxy, the slope S eoI (kpc" 1 km 2 s -2 ) of the Hi RC for the end of 
the EOR was compared among the 13 select galaxies. The Soor were found by 
fitting a least squares straight line to the outermost two or more data points of 
the RC with a 0.97 or higher correlation coefficient. The same procedure that 
was used to evaluate the i> 2 rmax - L relation was used to evaluate the S eoT - 
\K • a Q \ relation. The result is 



10 3 kpc _1 km 2 s" 



-7 B eor ± 15% , 



\K*d Q \ 

10 3 kpc _1 km 2 s" 2 



(25) 



where C* cmin = 0.96, K kl = 0.64±0.09, B kl = 1.74±0.06, and I soor = -4.8±0.1 
at la. Tables H] lists the k\ values for the 13 select galaxies. 

As seen in Fig. [2] the measured points in the OR appear nearly linear rather 
than a smoothly varying curve in most of the select galaxies. Therefore, the 
linear approximation of the RC slope in the OR was considered more appropriate 
than a smoothly varying function. Define S OY (km 2 s -2 kpc -1 ) as the slope from 
the -Rrrmax data point to the beginning of the S eoT of the Hi RC. For the 15 
select galaxies excluding NGC 2841, following the same procedure as for finding 
Eq. (Ull) yields 



^ = - K B h 

10 3 km 2 s- 2 kp C - 1 " h 11 lQSergs- 1 



lor ± 15% , 



(26) 



where the integer values are listed in Tabled! The C cm i n , K\ t , and values 
are listed in Tableland I or = 1.6 ± 0.1. 
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3.3 Second approximation 

Equation ([5]) indicates the parameters of galaxies depend on L and \K • a \. 
The first approximation selected either L or \K • <z | depending upon which is 
predominant. The second approximation considered the less predominant term 
to be a systematic correction term. Posit the residual between the measured 
parameter value and the calculated first approximation parameter value is pro- 
portional to the less predominant term. Appling the same procedure as for 
finding Eq. (fT2|) yields: 



^ rmax = K ai B a a ] + (-iy°K a2 B^ |A ^* Qo| 9 ± 10% , (27) 

10 3 km 2 s- 2 1 ai lOSergs- 1 1 ; aa a2 10 3 kp^ 1 km 2 s~ 2 1 ; 

where K a2 = 0.15±0.02, B &2 = 2.00±0.03at la, forNGC 5448 (<5v r 2 rmax /f r 2 rmax = -0.04), 
and for NGC 3031 (Kmax/^rmax = -0-03); 

= K B bi — + (-iy>K h2 B b b 2 \ K ' s o\ — 2% (28) 

kpc bl bl 10 s ergs" 1 v ; 2 b2 10 3 kpc- 1 km 2 s- 2 

where K h . 2 = 0.08 ± 0.01 and B h . 2 = 1.75 ± 0.08 at la; 

Rrrm^ = _L \ K * °° I ± 3% (29 ) 

10 3 kpckm 2 s- 2 1 Cl lOSergs- 1 V ; ° 2 C2 10 3 kp^ 1 km 2 s~ 2 1 ; 

where K C2 = 0.63 ± 0.07 and B C2 = 1.98 ± 0.03 at la; 

^rmax/^rrmax = * _L B * \K •a \ 

10 3 kpc" 1 km 2 s- 2 1 dl lOSergs- 1 V ' da da 10 3 kpc" 1 km 2 s" 2 ' 

where if d2 = 0.024 ± 0.007 and B d2 = 2.2 ± 0.3 at la; 

?E -K B ei L i I 1Y°K B e2 l-^'^l i n% fai) 

10 3 kpc km s -2 10 8 erg s 1 10 3 kpc km s~ 2 

where K C2 = 0.021 ± 0.003 and B C2 = 1.72 ± 0.03 at la; 

w 2 ^™ ov „f, L , I if • a, 



cormax 



TT^TTT + (-1)"^^ ' -iA_ 2 ± 2% , (32) 



10 3 km 2 s- 2 1 fl 10 s ergs- 1 v ; 2 fa 10 3 kpc" 1 km 2 s" 2 

where if fa = 0.7±0.1 and B h = 1.9±0.1 at la, for NGC 5448 (<5w 2 ormax /w 2 ormax = -0.03), 
and for NGC 3031 (<5« 2 orma > 2 ormax = -0-02); 

RZZZEL = Kgl B« T + (-iy°K g2 B£ \ K ' ^l % ( } 

kpc Sl 81 10 8 ergs" 1 v ; 82 82 10 3 kpc- 1 km 2 s- 2 

where K g2 = 0.008 ± 0.002 and B g2 = 2.05 ± 0.09 at la; 

^eormax^eormax " -r^ht ^ r \ e, -rr nlin \^ * 



10 3 kpc 

where if h2 = 7 ± 2 and B h2 = 2.06 ± 0.08 at la; 



^ormax/^cormax = 4l _£ (.^-ijf W> ^'" o| ? ± 2% , (35) 

lQSkpc^WV 2 1 11 10 8 ergs- 1 V ' 12 12 10 3 kpc" 1 km 2 s~ 2 V ' 



where K l2 = 0.024 ± 0.002 and B l2 = 1.78 ± 0.04 at la; 



Asym 7 ax = K h E* \ K *ao\ (-l)«K hBf'—s-^ T ± 8% , (36) 

10 3 km 2 s- 2 Jl Jl 10 3 kp C - 1 km 2 s- 2 V ' 32 ^lQSergs- 1 ' v ; 
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Table 7: Second approximation integer values for the select galaxies. 



Galaxy 


Oi 2 


b 2 


C-2 


d 2 


e-2 


h 




h 2 


>2 


7 9 


&2 


h 


NGC 0224 


























NGC 0598 


























NGC 3031 


8 


1 


7 


6 


4 


4 


3 


G 


1 


4 


3 


11 


NGC 0300 


























NGC 2403 


3 


3 


4 


4 


G 


2 


4 


2 


5 


5 


5 


6 


NGC 5457 


4 


G 


7 


5 


8 


3 


7 


7 


4 


5 


3 


8 


NGC 4258 


4 


2 


1 


3 


G 


2 


4 


3 


2 


4 


2 


1 


NGC 0925 


5 


2 


2 


2 


1 


1 


3 


1 


1 


4 


G 


3 


NGC 2841 












5 


8 


8 


4 


2 


4 


1 


NGC 3198 


6 


6 


9 


1 


2 


5 


8 


G 


4 


3 


2 


7 


NGC 4414 


4 


3 


6 


4 


G 


2 


5 


2 


1 


1 


1 


G 


NGC 3319 


4 


4 


6 


2 


4 


3 


7 


3 


5 


2 


G 


7 


NGC 7331 


7 


3 


6 


G 


7 


3 


G 


5 


5 


2 


3 


3 


NGC 4535 


7 


6 


5 


4 


9 


3 


4 


2 


7 


4 


3 


G 


NGC 4321 


1 


1 


2 


1 


2 


2 


1 


3 


3 


3 


1 


5 


NGC 4548 


16 


12 


8 


16 


21 


13 


14 


13 


16 


4 


5 
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where K h = 0.053 ± 0.008 and B h = 2.4 ± 0.2 at lrr [for the select galaxies, 
o- c ((5A symax /A symax ) = -0.04], for NGC 5448 (SA symax /A symax = -0.02), and 
for NGC 3031 (SA symax /A 

symax — 

-0.04); 

= K kl B^ (-l) s »K k2 B^ — ± 2% , (37) 

10 3 kpc- 1 km 2 s-2 kl lOSkpc^Ws- 2 V ' 2 ^lO^rgs- 1 ' V ' 

where K ka = 0.061 ± 0.006 and B k2 = 1.80 ± 0.05 at la; and 

*h B\] -J—, + J or + (-lr^X* J^'^s , ± 2% , (38) 



lO^m's^kpc- 1 1 11 10 8 ergs- 1 ° r v ; 12 12 10 3 kpc" 1 km 2 s~ 2 

where K h = 0.0039 ± 0.0005 and B h = 1.96 ±0.05 at la. If the measured value 
of a parameter is larger than the first approximation calculated value, the sign 
s x = in Eqs. (|27 p - (|3"8"|) . If the measured value of a parameter is smaller than 
the first approximation calculated value, the sign s x = 1 in Eqs. (|2"T)) - (l3"8"l) . The 
second approximation integer values are listed in Tables [21 and [7] 

The calculated RCs of the above equations are plotted as solid lines in Fig. [2] 
for the select galaxies. 

Because barred galaxies were part of the sample, the monopole Source sug- 
gests the p is not sourced by matter as is gravity. The Source appears as a 
monopole at the beginning of the RR that is a few kpc from the center of the 
galaxy. 



4 Discussion and conclusion 

The SPM L first approximation correlation with RC parameters is an expansion 
of the URC concept. The SPM includes observations such as rapidly declining 
RCs and asymmetric RCs. 

The second approximation equations show little effect of neighboring galax- 
ies on the Hi RC except when orientation and closeness of luminous galaxies 
produce a higher \K • a\. Although the neighboring galaxy's effect is less sig- 
nificant on the RC, it does induce small perturbations in the orbits of particles 
that cause the observed asymmetry. Because the m s /m L is smaller for higher 
metallicity particles, the effect of neighboring galaxies on the Ha RC is less. 
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Therefore, the finding of cluster effects on Ha RCs is subject to the sample 
selected and may be detected only when \K • a\ is large 

A large |-fT«a| requires the host galaxy be oriented appropriately and requires 
large Vp. A large Vp is present near the cluster center or where the near 
galaxies are large and asymmetrically arranged around the host galaxy. The 
former condition is caused by the r\ < of the galaxies in the cluster center 
surrounded by e > of the ga laxies in th e clust er shell. Therefore, the SPM is 
consistent with the finding of I Dale et al.l ( 2001 ). The latter condition may be 
caused by a large, near neighbor galaxy. 

The SPM's (GnigMg — G s m s L)/m L R 2 term has the form of a subtractive 
change in Newtonian acceleration. The SPM suggests the m s change with radius 
causes the species of matter to change with radius that causes the appearance 
of a gravitational acceleration change with radius that is related to L. If the 
|i4T»a| term is small, the SPM effe ctive force model reduces to the MOND model 
as was suggested in iHodgei ( 2006 ) for elliptical galaxies. 

The deviation of the data of NGC 5448 from Eq. (fT2|) suggest a physical 
mechanism behind the quantized galaxy parameters. The clear departure from 
circular mo t ion a nd the significant mass transfer inward (i? ^ 0) found by 
Fathi et al. ( 2005f ) suggests this galaxy is in transition from one virial state to 



another. Further, the noted stellar and gas velocity difference decreases at larger 
radii. The better fitting of the t>e 0rmax - L and of the j4 S ymax - |-K»a | relations 
is the expected result. N GC 3031 shows strong, non-circular motion in the disk 
( Gottesman et all . |l966). This suggests the integer variation is caused by the 
accumulation of mass at potential barriers such as at i?A and i? r rmax- Continued 
nucleosynthesis and changing |i^»a | causes an occasional, catastrophic rupture 
of one or more of the potential barriers, R ^ 0, and, therefore, the transition 
of the galaxy from one integer classification to another. A smoothly varying 
transition from the RR to the OR for flat or declining RCs such as NGC 4321 
suggests mass is accumulating at a potential barrier at the end of the RR and 
is being depleted from the outer parts of the OR. 

The proposition that observed asymmetries are caused by neighboring galax- 
ies may also apply to bars and to the formation and maintenance of the rotation 
of mass around the galax y. 

ISteinmetz fc Navarrol ( 20021 ) found in a series of N-body/gas dynamical sim- 
ulations that included feedback: that feedback is a necessary component for 
morphology determination; that the main morphological component is regulated 
by the mode of gas accretion and intimately linked to discrete accretion events; 
that morphology is a transient phenomenon; and that the Hubble sequence re- 
flects the varied accretion histories of galaxies. If luminosity is proportional to 
the e, which directly causes the parameters of the RC, then there must exist a 
feedback mechanism controlling the parameters of the RC. 

Approximately 66% of the sample galaxies have a\ — 4 or a\ = 5 as seen in 
Fig. [3] and Table [3] If a\ = 4.5 is used in Eq. (Ti"2"|) , v^ lmax a L for a majority 
of sample galaxies. Only NGC 4258 of the select galaxies would appear to 
be an outlier, which may suggest the w rr max is larger than the measured point 
chosen herein. Further, the neighboring galaxy effect would fail to improve 
the if rmax oc L relation. The effect of the integer values is to broaden the 
applicability of the parameter - L relations and to establish relations wherein 
the neighboring galaxy effect improves the calculation. 

The SPM does not use the mass-to-light ratio to fit RCs. Indeed, Eqs. (|14| 
and (fT9| imply the effective mass-to-light ratio varies among galaxies as ob- 
served. 

The use of the Mb derived from the TF and the degree of correlation achieved 
herein suggests the mass term is totally determined by baryonic matter. The Vp 
field serves the same role in the equations as DM. Equation (J3j) with the L term 
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and Eq. ([7]) with the Ma term suggests a much larger spiral galaxy baryonic 
mass [M of Eq. ([3])] than traditional Newtonian mechanics calculates(i?v 2 /G). 
Therefore, the Vp field causes Newtonian kinematic measurements to consid- 
erably underestimate mass in a galaxy [M w Rv 2 /G + G s m s L/Gm L ). The dif- 
ference between the SPM and other RC models of the "missing mass problem" 
is that the added mass follows the luminous mass as suggested by the m s /m L 
in the equations. Therefore, the added mass is baryonic rather than DM, a 
changing gravitation constant, a changing acceleration, or other form of "fifth 
force" . The "missing mass" is non-luminous, baryonic matter. 

The purpose of the present investigation was to expand the SPM to be 
consistent with galaxy RC, asymmetry and EOR observations. The resulting 
model adds the force of a scalar potential of the sample galaxy and of neighboring 
galaxies to the Newtonian rotation velocity equation. The form of the equation 
for each parameter of each region is the same with differing constants. Integer 
values, rather than scale factors, determine the particular galaxy parameter 
values. Among the sample galaxies, the B band luminosity of a galaxy is related 
(1) to the maximum rotation velocity, radius, mass, and acceleration parameters 
at the end of the RR and at the end of the EOR and (2) to the slope of the 
RC in the RR and OR. The scalar potential effect of the neighboring galaxies 
is related to the slope of the Hi RC in the EOR, to the asymmetry in the 
OR, and to the residual of the measured and calculated values of the above 
mentioned luminosity dependent parameters. The RC in the OR may be rising, 
flat, or declining depending on the galaxy and its environment. The scalar 
potential field causes Newtonian mechanics to considerably underestimate the 
mass in galaxies. The Source of the scalar field appears as a monopole at 
distances of a few kpc from the center of spiral galaxies. The resulting RC 
equation is consistent with the formation, evolution, and long term maintenance 
of asymmetry observations. The RC observations are consistent with the SPM. 
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A The many lines in the parameter relations are 
not random. 

The many lines in the plot in Fig. [3] suggest the data points may be random. 
The null hypothesis tested was "the data points are indeed random points". 
This null hypothesis was tested by following the procedure used in discovering 
Eq. (p~2|) as follows: (1) A trial consisted of generating 15 sets of two random 
numbers between zero and one and subjecting the sets to the slope, correlation 
coefficient, and a e tests the galaxy sample passed. Call the independent variable 
(Xi) and the dependant variable (Yi), where i varies from one to 15. (2) The 
equation tested was 

Y calci = K x B tl X t , (39) 
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where K\ is the proportionality constant, B is the exponential base determined 
by the relation to be tested, and U is the integer classification for the i set. 
(3) The Ki value for each trial was calculated as the minimum value of Ki = 
Yi/(BXi) with an Xi > 0.4. (4) The ti value was calculated for each set of 
(Xi, Yi) values as 



where "ROUND" means to round the value in the braces to the nearest integer. 
(5) If any subset lacked a value, that subset was ignored. (6) If the correlation 
coefficient of any ti subset of (Xi, Yi) values, including the (X , Yq) — (0,0) point, 
was less than 0.90, the trial failed. (7) The F C aici was calculated according to 
Eq. (|3"9"|) . If the cr > 0.21 between the Yi and Fcaici for the 15 sets, the trail 
failed. (8) A count of the number N n of sets of (Xi, Yi) with Xf + Y^ < 1 and 
the number A/xhy of sets of (Xi, Yi) with Xi < Yi was made. (9) The trial was 
redone with another 15 sets of random numbers for 30,000 trials. 

For B — 2 the results were: (1) The N„/(15 x 30, 000) = 0.78528 « tt/4 and 
Axity/(15 x 30,000) = 0.49999. Therefore, the random number generator per- 
formed satisfactorily. (2) Of the 30,000 trials 22 passed (0.07%), the remainder 
of the trials failed. Therefore, the null hypothesis was rejected with a 99+% 
confidence. For B = 1.65 the confidence level of rejecting the null hypothesis 
was greater than 0.95. 
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